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Abstract: An interesting problem in statistical physics is the condensation of classical 
particles in droplets or clusters when the pair-interaction is given by a stable Lennard- 

0^ . Jones-type potential. We study two aspects of this problem. We start by deriving a 

large deviations principle for the cluster size distribution for any inverse temperature 
/3 G (0,00) and particle density p G (0, pcp) in the thermodynamic limit. Here pcp > 
is the close packing density. While in general the rate function is an abstract object, 
our second main result is the F-convergence of the rate function towards an explicit 
limiting rate function in the low-temperature dilute limit /3 — t- 00, /O 4, such that 
J^ ' — /3~^log/9 — )• z/ for some u G (0,oo). The limiting rate function and its minimisers 

appeared in recent work, where the temperature and the particle density were coupled 
with the particle number. In the de-coupled limit considered here, we prove that 
just one cluster size is dominant, depending on the parameter u. Under additional 
assumptions on the potential, the T-convergence along curves can be strengthened to 

t^^ ] uniform bounds, valid in a low-temperature, low-density rectangle. 
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1. Introduction 



We consider interacting A^-particle systems in a box A = [0, LY" C W^ with interaction energy 



Un{xi,...,xn) ■■-- 



E 

i<i<j<N 



v{\Xi-Xj\), 



(1.1) 



where v. [0,oo) — t- M U {oo} is a pair potential of Lennard- Jones type. That is, 

• it is large close to zero, inducing a repulsion that prevents the particles from clumping, 

• it has a nondegenerate negative part, inducing an attraction, i.e., particles try to assume a 
certain fixed distance to each other, 

• it vanishes at infinity, i.e., long-range effects are absent. 

Additionally, we always assume that v is stable and has compact support. We allow for the possibility 
that f = oo in some interval [0, rhc] to represent hard core interaction. See Assumption (V) in 
Section Fl.21 below for details. 




Figure 1. The pair potential v{r) = 1.5r ^^ — 5r ^ of Lennard- Jones type. 

A particle configuration x = (xi, . . . , xj^) in the box is randomly structured into a number of smaller 
subconfigurations, that is, well separated smaller groups, which we call clusters. One of our main 
questions is about the joint distribution of the cluster sizes, i.e., their cardinalities. Intuitively, if 
the box size is large in comparison to the particle number, then one expects many small clusters, 
and if it is small, then one expects few large ones. We will analyse this question much closer in the 
thermodynamic limit, that is, keeping /3 G (0, oo) fixed and taking 



A^ 



oo. 



L = L 



N 



OO, 



such that 



Td 



P, 



for some fixed particle density p G (0, oo), followed by the dilute low-temperature limit 

1 



/3 — )• oo, p iO 



such that 



/? 



log/j 



1.2 



(1.3) 



for some v G (0, oo). In this regime, the total entropy of the system is well approximated by the sum 
of the entropies of the clusters, and the excluded-volume effect between the clusters as well as the 
mixing entropy may be neglected. As a consequence, particles tend to favor one optimal cluster size, 
which depends on v and may be infinite. 

In recent work [CKMSIO] , the free energy was analysed in the coupled dilute low-temperature limit 



A^ 



oo, 



/3 = /3. 



N 



oo, 



L = L 



N 



oo 



such that 



1 N 



:i.4) 



N 
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Figure 2 . A schematic figure illustrating the cluster decomposition of a particle configuration 
and the induced graph structure. 



with some constant v G (0, oo). It was found that the limiting free energy is a piecewise linear, 
continuous function of v with at least one kink, i.e., non-differentiable point. Furthermore, there was 
a phenomenological discussion of the interplay between the limiting cluster distribution and the kinks 
in the limiting free energy, on base of a variational representation. See Section [1.31 for details. 

In the present paper, we go beyond [CKMSIO] by considering the physically relevant setting of a 
thermodynamic limit and by proving limit laws for the quantities of interest. That is, our two main 
purposes are 

(i) to derive, for fixed /9,p G (0,oo), a large deviations principle for the cluster size distribution 

in the thermodynamic limit in (jl.2p . and 
(ii) to derive afterwards limit laws (laws of large numbers) for the cluster size distribution in the 
low-temperature dilute limit in ()1.3p . 



In this way, we decouple the limit in (jl.4p into taking two separate limits, and we prove limit laws for 
the cluster sizes in this regime. 

The organisation of Section [1] is as follows. In Section 11.11 we introduce our model and define the 
thermodynamic set-up. Our main result concerning the large deviations principle for the cluster size 
distribution is formulated in Section 11.21 The low-temperature dilute limit is discussed in Sections 11.31 
and ll.4[ Adopting additional, stronger assumptions we give in Section [1.51 bounds that are uniform in 
the temperature for dilute systems. Finally we discuss in Section [1.61 some mathematical and physical 
problems related to our results. 

1.1. The model and its thermodynamic set-up. Here are our assumptions on the pair interaction 
potential that will be in force throughout the paper. 



Assumption (V). The function v. [0,oo) 



. U {oo} satisfies the following. 



(1) V is finite except possibly for a hard core: there is a rhc > such that v = oo on (0, rhc) and 
V < oo on (rhc,oo). 

(2) V is stable, that is, Un{x)/N is bounded from below in A^ G N and x G ^'^d.'^^ 

(3) The support of v is compact, more precisely, b := supsupp(f) is finite. 
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v(r) 




Figure 3. Two examples of pair interaction potentials satisfying assumption (V). 



(4) V has an attractive tail: there is a 5 > such that v{r) < for all r G (b — 5,b). 

(5) V is continuous in [rhc,oo). 



Assumption (V) differs from Assumption (V) in [CKMSIO] in two points: here we drop the require- 
ment f (rhc) = oo, and stability was there a consequence of some cumbersome additional assumption. 

We introduce the Gibbs measure induced by the energy defined in (ll.ip . For /3 G (0,oo), A^ G N 
and a box A C M , we define the probability measure FL^^ on A by the Lebesgue density 

pW (dec) = _^-l-^--^e-^^^('^) da;, x G A^, (1.5) 



where 



Za{P,N)N\ 



1 



^^'^■"'^-M ./,„■= 



-/3C/Ar(cc) 



da; 



is the canonical partition function at inverse temperature /3. 

We introduce the notions of connectedness and clusters. Fix R G (6, oo). Given x = (xi, . . . ,X]s[) G 
(M'^)^, we introduce on the set {xi, . . . ,xn} a graph structure by connecting two points if their 
distance is < R. In this way, the notion of i?-connectedness is naturally introduced, which we also 
call just connectedness. The connected components are also called clusters. A cluster of cardinality 
/c G N is called a k-cluster. By Afc(a;) we denote the number of /c-clusters in x, and by 



PkM^) 



|A| 



the k-cluster density, the number of A;-clusters per unit volume. We consider the cluster size distribu- 
tion 

Pa := (Pfc,A)feeN (1-6) 



as an Mtv/iai -valued random variable, where 



Mp-.-- 



(P, 



feJyfceN 



G[0,oo)^| J]A;pfe<p}, pG(0,oo). 



:i.7) 



fcgN 



On Mp we consider the topology of pointwise convergence, in which it is compact. Note that for each 
finite A*" and any box A C M'^, 

N 
^kpk^A{x) 



k=l 



N 



a;GA^. 



However, some mass of p/^ may be lost in the limit A — )■ oo to infinitely large clusters. The distribution 
of p\ under the Gibbs measure Pi^^ is the main object of our study. 
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Introduce the free energy per unit volume as 

N 



/A(/3,^):=-^logZA(/3,iV). 

It is known |R99j that the free energy per unit volume in the thermodynamic limit, 

/(/3,p):= lim /[o,i].(/3,f^). (1.8) 

N/Ld-^p 

exists in M for all p > when there is no hard core, i.e., if rhc = 0. When rhc > 0, there is a threshold 
Pep > 0, the close packing density, such that the limit exists and is finite for p G {0,pcp), and is oo for 
p > Pep- Since we are interested in dilute systems, i.e., small p, we will always assume that p G (0, Pcp)- 

1.2. Large deviations for cluster distribution under the Gibbs measure. Our first main result 
is a large deviations principle (LDP) for the cluster size distribution under the Gibbs measure. For 
the concept of large deviations principles see the monograph |DZ98j . 

Theorem 1.1 (Large deviation principle with convex rate function). Fix j3 G (0, oo) and p G (0, Pcp)- 
Then, in the thermodynamic limit N — t- oo, L — t- oo, N/L — t- p, the distribution of p^ under Pa ^ 
with A = [0, L] satisfies a large deviations principle on Mp^i, with speed |A| = L , where e > is 
such that N/L < p + e. The rate function Jp^p'- Mpj^^ — t- [0, oo] is convex, and its effective domain 
{J/3,pi-) < oo} is contained in Mp. For p sufficiently small, {J/3^p{-) < oo} is equal to Mp. 

The proof of Theorem 11.11 is in Section [2l Define f{(3,p,-): Mp — )• [0, oo] through the equality 

JpM=:P{fiP,p,p)-fiP,p)). (1.9) 

Then the LDP may be rewritten, formally, as 

■^1 ^ e-''^-(^)]l{pA(a:) ^p}dx^ exp(-/31Al/(/3, p, p)] 

Thus f{(3,p,p) may be considered as the free energy associated with the cluster size distribution p^, 
thought of as an order parameter. The identity inf Jpp = translates into 

/(/?,p)=irif/(/3,p,.). 

Mp 

In words: the (unconstrained) free energy is recovered as infimum of the constrained free energy as 
the order parameter is varied, a relation in the spirit of Landau theory. 

It is a general fact from large deviations theory that an LDP implies tightness. More specifically, the 
LDP of Theorem 11.11 implies a limit law for the cluster size distribution towards the set of minimisers 
of the rate function. This is even a law of large numbers if this set is a singleton. Hence, Theorem ll.il 
gives us control on the limiting behaviour of the cluster size distribution under the Gibbs measure 
in the thermodynamic limit. However, in the general context of Theorem 11.11 we cannot offer any 
formula for the rate function Jp^p. We have to restrict ourselves to the low-temperature dilute limit 
(|1.3p . In this setting we obtain explicit asymptotic formulae in Section 11.31 below, and this is our 
second main result. 

1.3. The dilute low-temperature limit of the rate function. In this section, we formulate and 
comment on our main result about the limiting behaviour of the LDP rate function Jg p introduced 
in Theorem 11.11 and of its minimisers in the dilute low-temperature limit in ()1.3p . This behaviour is 
explicitly identified in terms of the ground- state energy of Um-, 

En := inf Un(x), A^ G N. 
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It can be seen like in the proof of |CKMS1^ Lemma 1.1] using subadditivity that the hmit 

Coo := lim — - G (-00, 0) 

exists. It hes in the nature of the regime in (II. 3p that it is not the cluster size distribution pk that will 
converge towards an interesting limit (actually, these will vanish), but the term g^ = kpf^/p, which 
carries the interpretation of the probability that a given particle lies in a /c-cluster. Therefore, let 

Q:={(%)fceMG[0,ir| J]<7fc<l} 

fcGN 

and introduce, for ly G (0, 00), the map g^'- Q — )■ M defined by 

9u{{qk)k) ■=^1k-^ \- (1- ^Qkjecx,- (1.10) 

fcGN fc6N 

Our second main result is the following. 

Theorem 1.2 (F-convergence of the rate function). Let v £ (0,oo). In the limit /3 — t- 00, p — t- such 
that — /3~^ logp — T- v', the function 



Q^MU{oo}, (g;-),^ !/(/?, p,(£2fe)^g^) 



T -converges to g^ 



For the notion of F-convergence, see the monograph |dM93j . Theorem 11.21 is proved in Section [5. 1[ 
The physical intuition is the following: at low density, the particle system can be approximated by 
an ideal gas of clusters, see |II561 Chapter 5] or |S03j . 'Ideal' means that we neglect the 'excluded 
volume', i.e., the constraint that clusters have mutual distance > R. As can be seen from the proof of 
Lemma Is. H this means that the rate function f{f3,p, •) is well-approximated by the ideal free energy 

f^°^'(/3,p,(p,)fc) ■.= ^kpkf^\(3) + [p-^kp,)f^^{^) + i^p,(logpfe - 1). (1.11) 

km fceN ^ km 

Here f^^iP) and f^{l3) should be thought of as free energies per particle in clusters of size k (resp., in 
infinitely large clusters), see Section [3] for the precise definitions. The functional pgi, is obtained from 
/ by two simplifications, justified at low temperatures. 



• First, we approximate cluster internal free energies by their ground state energies. 

• Second, we split the entropic term as 

-^pi,{logpk-l)=Y,Pk^ + ^^Pk[log^- 

^ fceN fcGN ^ ^ fcGN ^ 

and keep only the first sum. Thus we keep the entropic contribution coming from the ways to 
place the clusters (their centers of gravity) in the box and discard the mixing entropy. 

In classical statistical physics, the approach we take here goes under the name of a geometric, or 
droplet, picture of condensation |H561 [S03] . This is closely related to the well-known contour picture 
of the Ising model and lattice gases |R99j . Lattice gas cluster sizes have been studied, for example, 
in |LP77j . continuous systems were investigated in |M751 IZ08J . The focus of these works was on 
parameter regions where only small clusters occur. Our declared goal, in contrast, is to derive bounds 
that cover both the small cluster and the large cluster regimes (in the notation introduced below, this 
means both v > v* and ly < z^*). 

Under additional assumptions on the pair potential, we can replace the somewhat abstract F- 
convergence result with more concrete uniform error bounds, see Theorem 11.81 
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The rate function gi, appeared in |CKMS10] in the description of the behaviour of the partition 
function Zj^j^ in the coupled dilute low-temperature limit in ()1.4p . More precisely, it was shown there 
that, in this limit, for any v 6 (0, co). 

It was phenemonologically discussed, but it was not given mathematical substance to, the conjecture 
that the random variable q^^ = {kpk^AM/p)keN under Pa^\ with Ajy = [0,Ln]'^ satisfies an LDP 
with speed Nf3i\i and rate function given by giy{-) — ^i{^)- This would be in line with Theorem 1 1 . 1 1 and 
Theorem 11.21 and we do believe that this is indeed true, but we make no attempt to prove this. 

1.4. Limit laws in the dilute low-temperature limit. The minimiser(s) of the rate function 
/(/3, p, •) are of high interest, since they describe the limiting behaviour of the cluster size distribution 
under the Gibbs measure. It is a general fact from the theory of F-limits that F-convergence implies 
the convergence of minima over compact subsets and the minimiser(s). For the limiting rate function 
Qy, the global minimiser has been identified in |CKMS10] . The minimum is 

aiv) = migu = inf — — , (1-12) 

and the minimisers are given as follows. 

Lemma 1.3 (Minimizers oi g^). The number v* := mi]\f(zfq{E]\f — Neoo) is strictly positive. The map 
V I—)- iiiv) is continuous, piecewise affine and concave. Let M C (0, oo) he the set of points where fi{-) 
changes its slope. Then J\f is bounded, and /u(z^) = —i' for v > max A/" and \i{y^ = eoo for u < u* . 
Furthermore, 

(1) v* £ M C [i^*,oo), and Af is at most countable with u* as only possible accumulation point. 

(2) For V > V* , we have [x{y^ < Coo and every minimiser {qk)k of g^ satisfies Ylk&ilk — 1- ^/ 
V ^ A/", then g^ has the unique minimiser g'"' = {q)^ )k with q^ = 5]^Uy\ with k{u) the unique 
minimiser ofk>-^ {Ek — i')/k over'H. The map v i— t- k{y^ is constant between subsequent points 
in M. 

(3) For V < V* , we have ^{v) = Coo and the unique minimiser of gy is the constant zero sequence 
{(lk)k&n with qk = for any k. 

This is essentially [CKMSlOl Theorem 1.5], the proof is found in the appendix. If, as in [CKMSIO] . 
the point oo is added to the state space N of the measures in Q, then the minimisers of gi, are 
concentrated on N for u > i^* and on {oo} for z^ < z/*; it was left open in [CKMSlO] whether or not 
the latter regime is non-void. 

The set M is infinite if and only if {Ek — fceoo)fcGN has no minimiser. In dimensions d > 2, it is 
expected (and shown in some cases in, see [R81tlYFS09j ) that Ek — kcoo > est. k^'^'"^ — )• oo, ensuring 
that AA is a finite set. 

Now we can draw a conclusion from Theorem 11.21 about the limiting behaviour of the minimisers 
of the rate function in the dilute low-temperature limit. The following assertions are well-known 
consequences from the T-convergence of Theorem 11.21 see |dM93[ Theorem 7.4 and Corollary 7.24]. 

Corollary 1.4. In the situation of Theorem \1.S\ 

(1) the free energy per particle converges to fiif): 

^fif3,p)^^{u) 
P 
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(2) if IJ,{-) is differentiable at u (that is, for v G (0,oo) \N), any minimiser p'-'^'P'> 
Jpp converges to the minimiser of g^-' 

I.JP.P) 



{Pf^ of 



A; G N. 



, density g 










/ / -=^ " " " 


Phase IV 
(infinite 


/Phase III/ / 


configu- 


/ / / 


rations) / 


/phase II / 


^ 


^^^.^^ Phase I (single points) 

^fc^ 



temperature T 



Figure 4. A diagram illustrating the expected relationship of the slope condition 
— Tlogp = — /3^^ log/5 — 7- V and the minimisers of the rate function in the dilute low- 
temperature limit. 

Another important consequence of Theorem 11.21 together with the LDP of Theorem 11.11 is a kind 
of law of large numbers for the cluster size distribution p^^ in the thermodynamic limit, followed by 
the low-temperature dilute limit. A convenient formulation is in terms of the vector q^ = (gfc,A)fceN 
with g^^A = kpk^h/ P-, the frequency of particles in /c-clusters, if |A| = N/ p. 

Corollary 1.5. For any v G (0, oo) \M , any J^ G N and any e > 0, if (3 is sufficiently large, p 
sufficiently small and — 4logp is sufficiently close to v, then, for boxes A]\f with volume N/p, 



1™ JPjSivll^fcH-Aiv 



and 



N^oo 



hm P^ 



1| > e 



K 



Z] <iKhN > ^ 



if u > V* 



if V < V*. 



(1.13) 
(1.14) 



fc=i 



Proof. We prove ()1.13p and p.l4p simultaneously. Consider the set 



A 



p 



peMp-. 



peM,:Ef=i^>e 



e > for V > V* 
> for V < V* 



Then the F-convergence of Theorem 11.21 implies |dM93l Theorem 7.4] that 

liminf-irif/(/3,/>, •) > -mig^, 
I3,p p A A 

where liminf^^p refers to the limit in Theorem ll.2[ Furthermore, it is easy to see from Lemma [1.3l that 
S = inf ^ gi, — inf gi, is positive. We pick now (3 so large and p so small and — /3~^ log p so close to i' that 
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i inf A /(/3, p, ■) - inf A Qu > S/i and ^/(/3, p) - p-{y) < 5/4 (the latter is possible by Corollary [LiTl)!. 
Now the LDP of Theorem 11.11 yields that 

limsup-^logP)3^1^(pA^ G^) <-inf/;3,p = -/3[inf/(/3,p,-)-/(/3,p) 

< -f3p\mfg, - p{u) - I - ll = -/3p<5/2 < 0. 

LA -i -ij 

Hence, liniAr^oo 1P« ^ (PAn ^ ^) = 0. Noting that this probability is identical to the two probabilities 
on the left of (jl.lSp and (J1.14p for our two choices of A, finishes the proof. D 

It may come as a surprise that, for most values of the parameter z^, the cluster size distribution is 
asymptotically concentrated on just one particular cluster size that depends only on v. This may be 
vaguely explained by the fact that the zero-temperature limit /3 — t- oo forces the system to become 
asymptotically 'frozen' in a state in which every cluster size assumes the globally optimal configuration 
size, which is unique for v G (i^*,oo) \M. Furthermore, note that Corollarv 11.51 does not give the 
existence of 'infinite large ' clusters (i.e., clusters whose size diverges with A^) for any value of /3 and 
p, not even for u < u* and — /3~^ log/O w v. 

1.5. Uniform bounds. Under some natural additional assumptions on the pair potential, the asser- 
tions of Theorem 11.21 can be strengthened, see Theorem 11.81 below. Indeed, we will assume that the 
ground states of the functional Un consist of well-separated particles, which are contained in a ball 
with volume of order A^, and we assume some more regularity of the interaction function v. Then 
we show that the F-convergence in Theorem 11.21 in the coupled limit in (jl.3p can be strengthened to 
estimates that are uniform in some low-temperature, low-density rectangle (/3, oo) x (O,/?). This leads 
to corresponding uniform estimates on |-/(/3,/9) — /i(z^)| and on minimisers. We now formulate this. 

Assumption 1.6 (Minimum interparticle distance. Holder continuity). 

(i) There is rmin ^ ^hc such that, for all iV G N, every minimiser {xi, . . . ,xn) G (M ) of the 

energy function Un has interparticle distance lower bounded as \xi — Xj\ > rmin; i ¥" J- 
(ii) The pair potential v is uniformly Holder continuous in [rmimoo). 

The existence of a uniform lower bound rmin for ground state interparticle distance is, of course, 
trivial when the potential has a hard core rhc > 0. A sufficient condition for the existence of rmin > 
for a potential without hard core is, for example, that v{r)/r'^ — ?• oo as r — )• 0, as can be shown 
along [T06l Lemma 2.2]. 

Assumption 1.7 (Maximum interparticle distance). There is a constant c > such that for all 
N £ N every minimiser (xi, . . . ,xn) G (M'^)^ of the energy function Un has interparticle distance 
upper bounded by \xi — Xj\ < cN^''^. 

This assumption looks deceptively simple; on physical grounds, we would expect that it is true 
for every reasonable potential. To the best of our knowledge, however, non-trivial rigorous results 
are available in dimension two only, for Radin's soft disk potential |R81j and for potentials satisfying 
conditions (HI) to (H3) from |YFS09| . These potentials satisfy Assumption 11.61 as well. 

Theorem 1.8. Suppose that in addition to Assumption (V) the pair potential also satisfies Assump- 
tions lLd\ and \1.1l\ Then there are Ji^P^C > such that for every {P,p) G [/3, oo) x (0,p], putting 
u := — /3~^ logp, the following holds. 



(1) Estimate on the rate function: 
1 



C 



-f{P,P:{'-f)km)-9v{{qk)k) <^log/3> {qk)km^Q- (1-15) 
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(2) Estimate on the free energy: 

^f{P,p)-^{u) <2| log/3. (1.16) 

(3) Minimizers: For any niinimizer p'^'^-''^ of f{l3, p,-), put q^^-''^ := {kp^.'' /p)k£fi- Then, ifu < v* , 

EV^2-i—i log/3. (1.17) 

E ^^^^l-^A^y^log/^, (1.18) 



If V > V* , then 



where 

A(z^) := infj ^''"^ k G N\M(zy)| - p{v) > 

is the gap above the minimum, and M(z^) C N is the set of minimisers of {{Ej. — i^)/fc)fcGN 
(thus M{v) = {k{v)} for v ^ M). 

Theorem 11.81 is proved in Section 15.21 One can see from the proof that one can choose p = (2a + 
2R)~'^. It follows in particular that the F-convergence and the two convergences from Corollarv ll.4l can 
be strengthened to convergence for just taking /3 — )• oo, uniformly in p G (0,p], with an error of order 
/3~^ log/3. This form of the error order term is an artefact of the assumption of Holder continuity; the 
constant C depends on the Holder parameter. 

Note that (J1.17p implies that, in the case ly < v* , for every K G N, the fraction of particles in 
clusters of size < K is bounded by 

k<K ^ k<K ^ 

This shows that, as /3 — t- oo, for some choices of i^ = Ki^ — t- oo, the fraction of particles in clusters 
of size < Ki3 vanishes, i.e., the average cluster size becomes very large. Note that the law of large 
numbers in (J1.14p in Corollary [L5] may, under Assumptions 11.6] and 1 1 . 7| be proved also with K replaced 
byK/3. 

1.6. Some remarks concerning related mathematical and physical problems. Our prob- 
lem is connected with continuum percolation problems for interacting particle systems, see the re- 
view [GHMOl] . In our setting of finite systems, the term 'percolation' should be replaced with 
'formation of unbounded components', i.e., clusters whose size diverges as the number of particles 
goes go infinity. The problem of percolation or non-percolation for continuous particle systems in an 
infinite- volume Gibbs state (that is, in a grand-canonical setting) is studied in |PY09j . They prove 
that, for sufficiently high chemical potential and sufficiently low temperature, percolation does occur. 
However, they do not give any information on the densities at which percolation occurs. This hinders 
the physical interpretation, since one cannot say whether the percolation is due to high density or 
strong attraction. In this light, our results are stronger and at the same time weaker: we do show 
that a transition from bounded to unbounded clusters happens at low density, but only in a limiting 
sense along low-temperature, low-density curves; there is no fixed temperature or density at which we 
prove the formation of unbounded clusters. 

In addition, our work has an interesting relationship to quantum Coulomb systems. In the simplest 
case, a gas of protons and electrons, we may ask whether we observe a fully ionized gas, where protons 
and electrons stay for themselves, or a gas of neutral molecules, with protons and electrons paired up 
together. Rigorous mathematical results were given by |F85j . see also |CLY89j . in the Saha regime, also 
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called atomic or molecular limit: when the temperature goes to at fixed, negative enough chemical 
potential, the Coulomb gas behaves like an ideal gas of different types of molecules or particles. The 
chemical composition is determined by the chemical potential. 

Our results adapt this quantum Coulomb system picture to a classical setting. From this point of 
view, the key novelty is that we work in the canonical rather than the grand-canonical ensemble; this 
allows us to extend results to the region where formation of large clusters occurs. 



The remainder of this paper is organised as follows. In Section [2] we prove the LDP of Theorem ll.il 
in Section [3] we compare the rate function with an explicit ideal rate function, and in Section [J] we 
compare temperature-depending quantities with the ground states. Finally, the proofs of Theorems ll.2l 
and 11.81 are given in Section [5l 

2. Proof of the LDP 

In this section, we prove Theorem ll.il We fix /3 S (0, oo) and p G (0,/3cp) throughout this section. In 
Section 12.1! we explain our strategy and formulate the main steps, and in Sections I2.2H2.4I we prove 
these steps. The proof of Theorem 11.1! is finished in Section 12.51 



2.1. Strategy. The main idea is to derive first a large deviations principle for the distribution of 
{pk,k)k=\,...,j for fixed j G N, that is, for the projection of p^ on the first j components, and apply 
the Dawson-Gartner theorem for the transition to the projective limit as j — )• oo. From the proof 
of the principle for the projection, we isolate an important step, see Proposition 12.11 using standard 
subadditivity arguments, we prove the existence of thermodynamic limit for constrained free energy, 
the constraint referring to cluster size concentrations of size < j. The principle for the projection of 
p^ appears in Proposition 12.21 

Given N,Ni, . . . ,Nj £ No define the constrained partition function with fixed cluster numbers of 
size < j, 

1 r ^ 

ZAif3,N,Ni,...,Nj):=— e'^^^^") J] l{iVfc(^) = iVj da;. (2.19) 



fc=i 



Note that Za(/3, A^, A^i, . . . , iV,-) = if Ei=i ^Nk > N. 

In the following we shall often be interested in the interior or boundary of subsets A C [0, oo)-'^"^ 
for some j G N. Unless explicitly stated otherwise, Int A and dA refer to the interior and boundary of 
A considered as a subset of M-'"'^^. In particular, if G A, then is automatically a boundary point. 

We denote by dom/i = {x: h{x) < oo} = {/i(-) < oo} the effective domain of an (— oo, oo]- valued 
function h. 

Proposition 2.1. Fix j G N. Then there is a function fj{l3, •) : [0, oo)-'"''^ — t- M U {oo} such that 



fj{/3,-) is convex and lower semi- continuous; 

its effective domain has non-empty interior Aj := Intj^j+i dom/j(/3, •) and fjiP, •) is continuous 

in Aj,- 

its effective domain is contained in 

j 
dom fj{(3,-) C Aj C ^{p, pi, . . . , pj) G [0,oo)-'+^ pG [0,pcp], "^kpk < p|, 

k=l 
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and, moreover, if \An\, N, N['^ , . . . , A^- — )■ oo in such a way that 






then 



If(.P,Pi,---,Pj) G ^j, 

lim -^logZA^(/3,iV,iV{'^\...,iVf') = -/3/,-(/3,p,/5i,...,p,-)- (2-21) 

and i/ie /imii is finite. 

If {p, pi, . . . , pj) E dAj (boundary of Aj), then 

limsup-^ZA^(/3,iV,iVr, . . . , A^f ) < -(3fj{f3, p, p^, . . . , pj) G MU {-oo}. (2.22) 

If {p, pi, . . . , Pj) £ Aj , then (|2.2ip holds true and the limit is —(3fj{(3, p, pi, . . . , pj) = — cxd. 



This proposition is proved in Section | 

The set Ai is related to close-packing situations. For example, when j = 1 and the density p is 
higher than 1/\B{0,R)\ (where we recall that R is the parameter in our notion of connectedness), it 
is impossible to have a gas formed only of 1-clusters and we have fi{/3,p,p) = oo. 

Analogously to (|1.9p . let 

Il3,p,j{pl, ■■■, pj) ■■= Pi^fjW, P,Pl,---, pj) - fiP, P)) ■ 

We will prove in Section 12.41 the following. 

Proposition 2.2 (LDP for projection of Pa)- ^'^^ J ^ ^- Then, in the thermodynamic limit N — )■ oo, 
L —7- oo, N/L —7- p, the distribution of {pi^A, ■ ■ ■ j Pj,a) under the Gibbs measure 1P« a ^^'^ ^ = [0) ^] 
satisfies a large deviations principle with scale |A| and rate function I/i^pj- Moreover, the rate function 
is good and convex. 

Recall that a rate function is called good if its level sets are compact. In this case, it is in particular 
lower semicontinuous. The large deviations principle means that, for any open set O C [0, oo)-' and 
any closed set C C [0, oo)-'', with A = [0, L]'^, 

.rr^'^'d. ^logP<,"i((pi,A,...,p,,A)eC)) > -infJ^,,,,, (2.23) 

limsup ^logP;j~j((pi,A,...,PM)ec) < -inf//3,pj. (2.24) 

We refer to (j2.23|) as to the lower bound for open sets and to (j2.24p as to the upper bound for closed 
sets. 

2.2. Proof of Proposition 12.11 — subadditivity arguments. In this section we prove Proposi- 
tion 12.11 For the remainder of this section, we fix j G N. 

The crucial point is the following supermultiplicativity of partition functions, which translates into 
subadditivity of free energies: Let N' , N" G N. Let A', A" be two disjoints measurable sets which have 
mutual distance larger than the potential range b, and A large enough to contain the union of the two. 
Then 

Za{(3,N' + N") > Zj,,^j^„il3,N' + N") > Z^,{p,N')Z^„{p,N"). (2.25) 
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This standard trick leads to a proof of the existence of the thermodynamic hmit by subadditivity 
methods |R99j (where subadditivity is apphed to the microcanonical ensemble instead of canonical, 
but the method is the same). 

The starting point of our proof is the observation that a similar inequality holds for constrained 
partition functions Z\{f3,N,Ni, . . . ,Nj) provided A' and A" have mutual distance > R, where we 
recall that R G (6, oo) was picked arbitrarily. Therefore we can prove existence of the constrained free 
energy by adapting the standard methods. Let us recall, roughly, the standard strategy of proof: 

(1) As a first step, one proves existence of limits of —am log Z\{/3, N,Ni, . . . , Nj) along special 
sequences of cubes - roughly, the sequence is defined in an iterative way by doubling the cube's 
side length and adding a 'security margin', and multiplying particle numbers by 2'^. This uses 
subadditivity and yields a densely defined, convex function rj. 

(2) Then one shows that the function rj is locally bounded in some region of non-empty interior, 
and therefore can be extended to a continuous function / in some non-empty open set A. 

(3) At last, one proves the convergence of —-ahn log Za{P, N,Ni, . . . , Nj) to / along general cubes. 

Our proof follows these steps, with some complications. Notably, an extra argument is required 
in Step (2) (see Lemma 12.61 below). Moreover, we make the choice ~ convenient in view of the large 
deviations framework - to assign values to the free energy not only in A and outside A (where / is 
oo) but also in dA by requiring global lower semi-continuity and convexity. 

2.2.1. Convergence along special sequences. Let R' > R and Lq > be fixed, and define (L*)neNo 
recursively by L^^^ = 2L; + R' . Exphcitly, L; = -R' + T^{Ll + R'). Let A; = [0,^;]"^. Thus A*^^ 
can be considered as the union of T^ copies of A„ with a corridor of width R' between them. Let 

Vj := {p = {p,pi,. . .,pj) G [0,oo)^+i I p > 0, 3g G No : 2^^(LS + RYp G N^q""^}. 

Lemma 2.3 (Introduction of rjj{f3, •)). Let (p, pi, . . . , pj) G "Dj and put for n G N 

ATM := 2"^(L* + i?') V iV^ - 2"''(^o + ^') V (fc = 1, . . . ,i). (2.26) 

The following limit exists in M U {oo} and is equal to an infimuni: 



iyif3,p,pu...,p,):=-lim^-^logZ^,^if3,N^"\N[-\...,N^^-^ 



n— >oo 



inf f- -^ log Za. (/?, iV'"' , iVl"' ,..., iV'"0 
neNV B\A*\ "^ ' 1 ' ' J / 



(2.27) 



This limit is finite as soon as Za* (/3, N''"\ N^ , . . . , N^ ) > for some n G N. In particular, 

j 
{r]jiP,-)<(^}c[ip,pi,...,Pj)£Vj: ^kpk<p<Pcp}. (2.28) 

k=l 

Proof We can place 2'^ shifted copies of A* in A*_,_^ in such a way that the copies have distance > R' 
to each other. Hence we have 

Za.^^ (/3, N^-+'\N[-+'\ ..., Arj"+^') > (Za. {(3, N^-\N[-\ ..., iVJ"'))' . 

Abbreviating 

1 9^1 A* I 



— log Za* (/3, A^("' ,N["\..., N"^) and 1 + e„ : 



|A*| O ^^-nV"' ' i ' ^ J / ' " I A* I 
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this is just the inequahty Un+i < (1 + £n)un- Our goal is to show that hm„_s>oo Un exists and is equal 
to n := inf„gN Un- Remark that 

_2d\K\ _f 2r^+HL* + R')-2R' Y_ 

which yields ^^^ \£n\ < oo. The case u = — oo is excluded by exploiting the stability of the energy: 
for some C € (0,oo), we have 

ZA,(/3,iV('^),<',...,iVf' ) <^A,(/3,iV("') < ^e-^V") |A:r^"' < e^^'"\ 

and hence u > —Cp. 

If u = oo, then ti„ = oo for all n and in particular n„ — t- oo = u- Consider now the case u G M. For 

5 > 0, let gr G N such that Ug < i + 6 and 1 - 6 < 111=^(1 +ek) <l + 6 for ah n > q. Then for n > q, 

n-l 
U<Un <UqJ\{l+ek) < {u + 5){l+5). 

k=q 

Letting first n — )• oo and then 5 — )• we conclude that Un — >• u- The additional assertion is clear from 
the proof and from the fact that, for p > pcpi we have oo = /(/3,p) = — 4lim„_^oo tttt log^A* (/?, A^'"^). 

D 

2.2.2. Properties of the limit function r]j{/3,-). The next lemma essentially states that r]j{(3,-) is a 
convex function. The precise formulation needs some care since the domain T)j of this function is not 
closed under taking arbitrary convex combinations. 

Lemma 2.4. Let p,p' G T)j. Let t G (0, 1) he a dyadic fraction, i.e., of the form t = p/2'^ for some 
p,q e No- Then tp+ {1 - t)p' G Vj and 

rjj{/3,tp+{l - t)p') < tr]j{P,p) + (1 - t)r?,(/3,p')- (2-29) 

Proof. Consider the cubes A* defined as above. A*_,_]^ is the union of two sets of 2'^~^ copies of A* 
plus some margin space. We first consider i = i. We can lower bound 

ZA;^^(/3,2("+i)'^(LS + i2')'(P + P')/2) 

> (zA.(/3,2"'^(L5 + i?')V))''"(^Aj/3,2"'^(L^ + iZ')V))'' '• 

We divide by |A*_,_]^| and pass to the limit, this gives Eq. (j2.29p for the case t = ^. The general case 
is obtained by iterating the inequality. D 

The following is a technical preparation for the proof of the local boundedness of r/j (/?,•) in 
Lemma 12.61 and will also be used later. We define a cluster partition function with volume constraint: 
for a, /3 > 0, /c G N, let 

Zl''^{/3) :=-!-: f e-'^^("^'"2'-'"'=)l{{xi,X2,...,Xfc} connected} dxi---dxfc. (2.30) 

Lemma 2.5. Let 6 G (0, [R — rhc]/3). There is a C{6) G M such that for all k £ N and a^ > 

6 + e/'^{rh, + 26), 

4zt'''\P) > \B{Q,5/2teM-PC{5)k). (2.31) 
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Proof. The cube [0, Uk]'^ is large enough so that, for some h G (rhc + 26,R — 6) and some 9 G M'^, the 
cubic lattice [0, a^J'^n [9 + (/iZ)'^) contains at least k points all having distance > 6/2 to the boundary 
of the box. By placing particles in the lattice, we obtain an (R — 5)-connected reference configuration 
(xi, . . . , Xk) S ([0, OkY)'' with the following properties: 

• All points have distance > 5/2 to the boundary of [0, a^]'^. 

• Distinct points Xi,Xj have distance > r^c + 6 to each other. 

We can lower bound Z^ '"''' {(3) by integrating only over those configurations with exactly one particle 
per ball B{xi,6/2). Such a configuration is always i?-connected. Moreover the energy of such a 
configuration can be upper bounded by C{6)k with 

C{6):= Y^ sup |7;(s|£|)| <oo, 

and Eq. dOT]) follows. D 



Lemma 2.6 ({??j(/3, •) < oo} has non-empty interior). For p G (0, oo), let 

j 
^jip) ■= {{p,Pi,---^Pj) e (0,oo) X [0,00^ p<p, Ykpk< py 

k=l 

Let 5 G (0, {R - rhc)/3) and C{5) he as in Lemma \KR Fix p{5) := (rhc + R + 25)~''- . Then for all 
p G Aj{;p{5)) n Vj, we have r]j{(3, p) < C{6) - /3~^ log \B{0, 6/2)\ < oo. In particular, 

A,{p{6))nV,c{rjjW,-)<oo}. 

Proof. We first give an appropriate lower bound for the constrained partition function for the two 
extreme cases when (1) all clusters have the same size k G {1, . . . ,j}, and (2) all clusters are larger 
than j. Afterwards, we use the convexity of ?7j(/3, •) (see Lemma l2.4p to handle all other cases. 

Thus fix p = {p, pi, . . . , pj) G Vj n Aj{'p{6)). In the first case, let k G {1, . . . , j} and p = p'*"' with 
p(f) = Pk = p/k and pf^ = p^ = for i / A;. It follows that the A^'"', Af^"''s defined as in Eq. (|2:26D 
satisfy A^*"' = kN^"'' and A'^"' = for i / fc. Furthermore, let Ok > 5 + k^/^^ir^c + 26) such that 
p{ak + R)'^ < k. We are going to use the boxes A* defined above. In A*, we place cubes of side-length 
Ofc with mutual distance > R. As n — )• (X), the number of such boxes behaves like 



ir, 



\A*\ 



{ak + RY 



> 



(ok + Ry k 

Thus we can lower bound the partition function by requiring that each /c-cluster lies entirely in one 
of the above boxes, and there is at most one cluster in each such box. This gives 

ZA.(/3,Ar("),Arr\...,Arf ) > (^(^'i/^) (af Z,^^''^'=(/3))'^'"'^' > |B(0,5/2)r^"^xp(-/3Ar(")C(<^)), 

(2.32) 
where in the last step we used Lemma [2 .51 and estimated the counting term against one. Thus we find 

lim -^logZA. (/3, A^("', Ar{"', . . . , A^j"') > p(-PC{5) + log \BiO,6/2)\) . 

n—^oo I A* I " J- ./ \ y 

Thus, 

Vj{l3,p''') < p{c{5) - r'log\B{0,6/2)\y 

In the next step, we assume that p = p'*^^ with pj° = p^ = for all k = 1, . . . , j. Again, we define 
N^"^ and the NJ^"^ by (|2.26p . We now lower bound the constrained partition function by putting all 
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particles into one cluster. 

ZA.(/3,iV("',iV|"\...,Arj-') > |A:|z;;'(5(/3) for iV^"' > j + 1. 
Observe that a„ '■= L'^ satisfies the conditions from Lemma 12.51 thus we also have 

r/,(/3, p*"') < p{c{5) - r' log 1^(0, <5/2)|) • 

In the general case, let q^ := kp^/ p for A; G {1, . . . , j} and q^ := 1 — X^^.^^ Qk- Then qq, qi, . . . , Qj > 
are dyadic fractions and satisfy J2k=o1k = 1- Furthermore, p = '^i^oQkP'^''^- It follows from 

Lemma 12.41 that 

j 

r/,(/3,p)< j;(?fcr/,(/3,pW)<p(c(5)-r'log|B(0,V2)|). 

fc=0 

D 

2.2.3. Extension of rij{l3, •) to W~^^ . We now extend ?7j(/3, •) : Pj — )• M U {oo} to a convex, lower semi- 
continuous function /j(/?, •): M-'+i ^ MU {oo}. We follow the proof of [R99l Prop. 3.3.4, p. 45]. Let 
Tj be the closure of {r/j(/3, •) < oo}, and let Aj be the interior of Tj. Note that Tj C [0, oo)-'"'"^, as 
r/j(/3, •) = oo on M^'+i \ [0, 00^'+^ 

Lemma 2.7. (1) The interior Aj ofTj is non-empty. 

(2) The restriction of r]j{(3, •) to Vj n Aj has a unique continuous extension fj{l3, •) : Aj — )• M. 

(3) De/ine /j(/3, •) : M^+i ^ R U {00} by 

fjW,p) ifp^^j, 

+00 ifpeA'j, (2.33) 

hminf^,^, /,(/3,p') ifpedAj. 



f0,P) 
Then fj{f3,-) is convex and lower semi- continuous, ana 



p'eA^ 



fj{P,p) = lymf,{l3,p + t{p' - p)), p G 5A„ p' G A,. (2.34) 

(4) 

i 
{/,(/3, •)< 00} C A, c{(p,pi,..., /5j)G [0,00)^+1 I pG [0,pcp], J]fcpfe<p}. (2.35) 

fe=i 

Proof. (1) This follows from Lemma 12.61 

(2) For the existence and uniqueness of a continuous extension in Aj, follow \R99\ p. 45]. The key 
point is that in Aj, i]j{j3,-) is a locally uniformly bounded, densely defined, convex function in the 
sense of Lemma 12.41 

(3) Let us extend Jj{(5, ■) to W+^ with /^(z?, p) = 00 for p G MJ+^\Aj. Then /,(/3, •) is convex, but 
may fail to be lower semi-continuous. Furthermore, /j(/3, •) and /j(/3, •) can differ only on i9Aj. The 
lower semi-continuous hull of /j(/3, •) is 

dJj{P,p) := liminf fj{P,p'), p G W+\ 

see |HLOH Def. 1.2.4, p. 79]. This is a convex, lower semi-continuous function which coincides with 
fj{l3,p) in Aj jHLnil Prop. 1.2.6, p. 80]. It follows that cl/j(/3,p) coincides with /j(/3,-) in Aj. It 
is elementary to see that in the definition of cl/j(/3,-), the limit inferior can be restricted to those 
p' ^ p that are in Aj. In other words, cl/j(/3, •) and /j(/3, •) coincide throughout M-^^^. This shows 
that /j (/?,•) is convex and lower semicontinuous. Eq. (|2.34p follows from [HLOlt Prop. 1.2.5]. 
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(4) The first inclusion follows from the definition of /j(/3, •), and the second from (I2.28p . D 

2.2.4. Limit behavior along general sequences. 

Lemma 2.8. Fix {p, pi, . . . , pj) G (0, oo)-^"*"-^. Let (A^-[^')7vgM5 • • • , (-^i^'')7VgM be No-valued sequences 
and {Aj\f)]\[^^ a sequence of cubes such that as N ^ oo, (j2.20p holds. Then, if {p, pi, . . . , pj) is in Aj, 

hm -^ log Za^(/3, N, iVf ', . . . , N'/') = -I3fj{l3, p, pi, . . . , pj) eR. (2.36) 

Proof. We proceed as in |R991 pp. 47]. We first prove the lower bound in (j2.36p . We will approximate 
{p, pi, . . . , Pj) with (p* , p\, . . . , p*) £ Vj satisfying p* > p and p\ < pi, ■ ■ ■ , p* < Pj- The idea is to pick 
the size parameter n = n{N) — )• oo of the special sequence of cubes A* z^-, introduced at the beginning 
of Section 12.2.11 in such a way that the cubes are small compared to Ajv- Hence, we can place a 
lot of them inside A^v at mutual distance > R. Afterwards, we distribute the particles and clusters 
inside a certain number of special cubes according to the distribution {p*,pl, . . . , p*) and place the 
few remaining particles somewhere else in A^v- 

Let {n{N))iyi=^ be an integer-valued sequence such that 

n{N) -^ oo and |A;(^)|V|AAr| -^ 0. 

We define Ni and NJ'j^'^ by (12.260 with n replaced by n{N) and p,pi,...,pj replaced by 

P*,pI, . . . , P*. Let ruN G Nq' and r^^' G {0, . . . , A^i"'"^" - 1} be such that 

This is possible because p > p* and therefore N > Ni for all sufficiently large N. For k G 

{l,...,j}, define rj!'^ by 

We claim that, for sufficiently large A^, the rj^ are non-negative integers. Indeed, this follows from 



*| A* I Pkl^MNV ~ a* 
P \^\i{N)\ P 



A7(iV) lA I J Ar("(-'V)) F -"-JV *ia* I *IA 



in combination with Pk ^ pI- > -^9%- Moreover, we can place m^ + r'^' copies of A* ,^n with mutual 
distance > R inside A^v- This is so because 

mN\KiN)\ ~ J^^N\ and r(-'|A;(^)| = 0(7V("(-»|A;(^)|) = 0{p*\KiN)?) = o{\An\). 

We lower bound the constrained partition function with parameters N, Nj^ \ . . . , A'^- by distributing 
first particles and clusters in the mj\f boxes following the distribution iV^"^'^''. This leaves r^'^' particles. 
Of those we distribute first k rj^ as clusters of size k, one per special cube, and then we distribute 
the remaining s*'^' particles into clusters of size i + 1 except maybe for one of size between j + 2 and 
2j + 1. Pretend for simplicity that they all have size j + 1. Then we get 

logZA,(/3, Ar,iv(-), . . . ,iv(-)) > m^logZA.^^^(/3,iVi"W), Ar(;(-» . . . ,iV3"») 

i+i 



k=l 



where L*,^^ denotes the side length of A*/^^. Using that Yl''k=i''^k'^ — ''''^' — A^i" ''^" = o(|AAr|), we 
get 

liminf-^logZA^(/3,7V,Ar{'^', . . . ,iVf ') > -/34/,(/3,p*,pt, . . . ,p*). 
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Now let {p* , pI, . . . , p*) — )• {p, pi, . . . , pj) and use the continuity of fj{/3, •) in Aj, to obtain 

lim inf -^ log Za^ (/3, N, N[''\ . . . , N^ ) > -/3/,- {P, p, pi, . . . ,pj). 

Now we prove the upper bound in (j2.36p . First of all, let us observe that the lower bound holds 
not only for sequences of cubes, but more generally for sequences of domains A^ that converge to 
infinity in the Fisher sense, as can be shown along the lines of our proof and |R99j . We shall need 
the statement not for general Fisher domains but only for A^ defined below, which is an L-shaped 
domain that is a difference of two cubes. 

Now fix C G (0, ^). For iV G N, let n{N) G N be so large thatAL^N contains An and satisfies 



0<C<-|^<^, nGN. 

I n(Nj' 



Let A^ be the set of points in A\^. having distance > R' to A^. Then (|AAr| + |A'^|)/|A*.jyJ -)■ 1. 
Let p* = {p*,pl,...,p*)€ Aj n Vj such that p* > 0. Define 7Vi"<^» and A^'J^" as in Eq. (I236D with 
n replaced by n{N) and p, pi, . . . , pj replaced hy p* , pi, ..., p*. Then 

^A;(,,(/3,A^r"",A^i:r\---,A^i?"") 

(2.37) 
Assume for simplicity that |AAr|/|A*j-^J — )• a G (0, 1/2] (otherwise go to suitable subsequences). Then 



ATJ"'^" - N P*\K{N)\ - P\^M p*-pa 



n 
--■■ P ■ 



Define /?'/,..., p"- in an analogous way and put p" = (p", p'(, . . . , p"-). Thus p* = ap + {\ — a)p" and 

\p -p\ = (i-a) \p- p \<2\p- p\, 

with I • I the Euclidean norm. Let e > such that Bir{p) C Aj. Now additionally assume that 
p* G B^i2{p)- Thus p" G Aj. In Eq. (12.370 . we take logarithms, divide by |A*.jyJ and pass to the 
limit A^ — 7- oo, which gives 

-/3/i(/3,p*)>alimsup-^logZA^(/3,Ar,Arf\...,ArfV(l-«)/3/j(/3,p")- 
To conclude we let p* ^ p (hence p" — > p) and use a > and the continuity of /j(/3, •) at p. D 

Lemma 2.9. Assume the situation of Lemma \2.8[ If p = {p, pi, . . . , pj) is in A • or in dAj, then 

hmsup -^ log Za^(/3, A^, Af{^\ . . . , A^f ) < -/3/j(/3, P). 

(Recall that fj{f3,p) = oo in the first case.) 

Proof. We proceed as in |R991 Prop. 3.3.8, p. 48]. One can show that there is an a G (0, 1/2] such 
that for p* G T>j satisfying p^ > whenever p^ > 0, and p" G Aj satisfying 

p* = ap+(l-Q)p", (2.38) 

it holds that 

- /37?j(/3, p*) > ohm sup -r^ logZ^J(3, N, A^f \ . . . , N^) - (1 - a)/3/j(/3, p")- (2-39) 

The proof of this is similar to the proof of the upper bound in Lemma [ 



Qlimsup-— logZA^(/3,iV,iVr,...,iVf')<limsup -/3/,(/3,p*(e)) + (l-a)/3/,(/?,p"(e)) 
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a) Consider the case p (^ Aj. For p" G Aj, we define p* by (|2.38|) . By choosing p" close enough to 
dAj , we can ensure that p* G Vj n A • . Thus we conclude from (j2.39p that 

lim sup -^ log Za^ (/3, N, iV{~' , . . . , N^^''^) = -oo. 

b) If p G dAj, let p'(e) G A^ n ^^(p) be such that fj{p,p'{e)) -^ fj{l3,p) as e i 0. By [HLOTl 
Lemma 2.1.6, p. 35], the half-open line segment (p, p'(e)] is contained in Aj. Since T>j is dense and 
because of the continuity of /j(/3, •) at p'{s), we can find p"{e) G Aj n Bi,{p) such that 

• p*(e), defined by dOSjl with p" replaced by p"(e), is in Aj n Pj n -Be(p). 
. |/i(/3,p'(e)) - /i(/3,p"(e))| < e, so that /,(/3,p"(e)) ^ /,-(/3,p) as e ^ 0. 

It follows from Eq. dOO]) that 
1 

TV-i-oo I^ATl - - ■' e_^o 

<-a/3/,(/3,p). 

D 

Proof of Prop. \2.1\ This is now straightforward from the previous lemmas. D 

2.3. The p-sections of Aj. We already saw that the set {/j(/3, •) < oo} has non-empty interior Aj. 
In view of the large deviations principle we are interested in properties of the map (pi, . . . , pj) i— t- 
fj{P,p,pi,...,Pj) at fixed /3 and p. This means that we look at the restriction of fj{/3,-) to the 
hyperplane of constant density p. 

Now, this restricted map inherits the convexity and lower semi-continuity from /,(/?, •). The question 
whether the set where it is finite has non-empty interior is, however more subtle. Closely related is 
the question whether Aj has non-empty intersection with the hyperplane of constant density p. 

To this aim consider the p-section of Aj, 

Cj(p):={(pi,...,pj)G(0,ooy I (p,pi,...,pj)GA,}. (2.40) 

Put differently, {p} x Cj{p) is the intersection of Aj with the hyperplane of constant density p. The 
hyperplane always cuts through the interior of Aj, i.e., cannot be tangent to Aj: 

Lemma 2.10. For any p G (0, pcp); t^^ set Cj{p) is non-empty, convex and open. Moreover, 



C,{p) = {(pi, ...,pj)e [0,oo)^ I (p,pi, . . . ,p,) G Aj}. (2.41) 

This last equation says that it does not matter whether we take first the p-section and then close 
the set, or if we close first and then take the section. 

The essential ingredients of the proof of Lemma 12.101 are the convexity of /j(/3, •); Lemma 12.61 and 
the following. 

Lemma 2.11. Let p G (0,pcp). Then there is at least one point (pi,...,pj) G [0,00)-^ such that 
fjW,P,Pi,---,Pj) < oo- 

Proof Let iV/|Aiv| ^ p. Let {N['^\ ..., TVJ"^') be such that 

Za^ (/3, N, N[''\..., iVf >) = max Za^ {f3,N,Ni, . . . , N^). 

{Ni,...,Nj)mi 

According to the Hardy- Ramanuj an formula, the number of partitions of N is not larger than 
exp(0(ViV)). Thus we find 

^A^ (A iV) < exp(0(ViV))ZA^ (A AT, iV}'^' , . . . , ATf ). 



20 SABINE JANSEN, WOLFGANG KONIG, BERND METZGER 

Passing to a suitable subsequence, we may assume that NJ^'^^/\An\ — )• pk, k = l,...,j, for some 
(/9i, . . . , pj) £ [0,00)-'. The previous inequahty then yields 

-00 < -/3fi/3,p) < -f3fj{f3,p,pi,...,pj). 

n 

Proof of Lemma \2. 1 (A Let p G (0, pcp)- Let p' G {p,Pcp) and [p'^, . . . , p'-) G [0, cx))-' such that 
fj{j3,p') < 00, where p' = {p',p'i, ■ ■ ■ ,p'j)- Hence, p' G Aj. Let p(5) and A{'p{6)) be as in Lemma [2T6l 
Let C C [0, 00 )■'"'' ""^ be the cone with apex p' and base A{p{6)), i.e., the set of convex combinations of 
points in A{'p{5)) and p'. By convexity, C C Aj. Looking at the p-sections of C we find that Cj{p) is 
not empty. 

Convexity and openness of Cj{p) are inherited from Aj. 

Now we prove ()2.4ip . Let Hp = {p} x W C W^^ be the hyperplane of density p. By jHLOU 
Prop. 2.L10, p. 37], 

Aj n Hp = A]nH^. 



The left-hand side is identified as {p} x Cj{p) while the right-hand side is {p} x A with A the set from 
the right-hand side in Eq. (|2.4ip . D 

2.4. Proof of Proposition [2T2] LDP for the projection of p^. In this section, we prove the large 
deviations principle for (/9i^A; • • • )Pj,a) under the Gibbs measure, as formulated in Proposition 12. 2[ 
This is equivalent to showing the two bounds in (j2.23p and (j2.24p and the claimed properties of 
^l3,P,j- Observe that the distribution of {pi^A, ■ ■ ■ ,Pj,a) under the Gibbs measure is concentrated on 
the compact set Mp. Hence, the family of these distributions is in particular exponentially tight. 
Hence, it is enough to prove the upper bound in (I2.24p for compact sets. From this, in particular the 
compactness of the level sets of I(3,pj follows, but we will also give an independent proof. 

For the remainder of this section, we fix p G (0, Pcp)- 

2.4.1. Properties of Ijj^pj. Recall the function 1/3, pj'- [0, oo)-' — )• MU {00} from ()2.ip and the p-section 
Cj{p) of Aj from ()2.40p . Recall from Lemma 12.101 that Cj{p) is non-empty, open and convex. 

Lemma 2.12. (1) Ip^pj is convex, and its level sets are compact. 

(2) I/s^pj is finite in Cj{p) and infinite in the complement of the closure of Cj{p). 

(3) For every open set O C [0,oo)-', 

infl« .= /i^fc?nc,{p)^/3,p,i ifOr\Cj{p)^0, 
o ^'"'^ [00 i/OnC7j(p) = 0. ^ ■ ' 

Remark 2.13. Eq. (j2.42p will he needed in the proof of the lower bound for the large deviations 
principle. The convexity enters in a crucial way in Eq. (j2.42p . Lower semi- continuity alone would not 
suffice! - (3) proves that the open set Cj{p) is a L 13 ^pj- continuity set, see JDZ981 p. 5]. 

Proof. (1) Convexity and lower semi-continuity are immediate consequences of the properties for 
/j(/3, •), since the restriction of a convex, lower semi-continuous function to a hyperplane is also convex 
and lower semi-continuous. Thus the level sets of I/3,pj are closed. By Eq. (j2.35p . 

j 
{Il3,p,j < 00} C |(pi,...,pj) G [0,00)-' I '^kpk < pj. 

k=l 

It follows that the level sets are also bounded, hence compact. 

(2) If {pi,...,pj) is in Cj{p), then (p, pi, . . . , pj) G Aj by definition of Cj{p), and therefore 
/j(/3, p, pi, . . . , Pj) < 00. Hence, Ii3,p,j{pi, . . . , Pj) < 00. 
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If {pi, ... , pj) is in the complement of the closure of Cj{p), then by Eq. ()2.4ip . {p, pi, . . . , pj) is in 
the complement of the closure of Aj, from which I/3^pj{pi, . . . ,pj) = cx) follows. 

(3) If O and Cj{p) are disjoint, I/s^pj = +00 on O by (2). If the sets are not disjoint, we know that 

i",f ^/3,P,i = 'BL_ h,P,j < ^^ inf Ip^p^j , (2.43) 

and it remains to prove the opposite inequality. Thus let p = (pi, . . . , pj) G OridCj{p). Let p' £ Cj{p). 
By Eq. ^Sl, 

I/s,p,j{p) = l™^/3>p>i(^ + *(^' - ^))- 

Because O is open and by [HLOll Lemma 2.1.6, p. 35], for sufficiently small t, p+t{p' — p) € Or\Cj{p). 
Thus for some suitable to > 0, 

Ip,P,j{p) = lim//3,pj(p + t{p' -p))> inf I^,pj(p + t{p' - p)) > inf I^^p^^ 

t^yj i€(U, 10} U\\Cj{p) 

u 

2.4.2. T/ie two hounds in ^TM and ^a?2M . For A C [0,oo)-'', let 

3 

VN{j,A):=[{Ni,...,Nj)eni\ (iVi/|A;v|,...,iV7v/|Av|) gA Y^kN^ < n]. 

k=l 

We note that the probability of finding (pi,Ajvi • • • 1 PjAn) i^ ^^^ ^^^ ^ is a sum of constrained partition 
functions: 

^'p'',L{iPiA.,---,Pj,A.) ^^) = z, (B N) ^ Z^^{f^,N,Ni,...,Nj). 

^^^^' > {Nu-,Nm)£Vn(jA) 

Upper bound in (|2.24p for compact sets. Let K C [0,oo)-' be a compact set. Let {N^ , . . . ,N-) G 
Vn{J,K) maximize the constrained partition function over VnU^K), i.e., 

ZA^(/?,iV,iVf\...,iVf')= max Za^{(3,N,Ni, . . . ,N,). 

{Ni,...,Nj}eVN(j,-K) 

Then 

«^ ((P1.A., • • • , P.A.) e i^) < ^|;§y^A^(/3' ^^ <^ • • • ' ^^r )• 

Now, the cardinality of Vn{j,K) is smaller than the number of partitions of N, and therefore not 
larger than exp(0(\/iV)), which is e"'^^^). The sequence (A^{^'/|AAr|, . . . , A''!^'/|AAr|)7VGN takes values 
in the compact set K and therefore, going to a subsequence, we can assume that it converges to some 
{pi, . . . , Pj) G K. Applying Proposition 12. II we ffiid 

limsnp-^logZA^{P,N,N[''\...,Np)<-f3fj{f3,p,pi,...,p,)<-f3mffj{l3,p,-). 
This yields the upper bound in ()2.24p for K = C. 

Lower bound in ()2.23p for open sets. Let O C [0, oo)-' be an open set. Let {pi, . . . , pj) € O. We can 
choose {N[''\ ..., A^j^^) G Vn{J, O) so that A^^^VjAtvI -^ Pk, k = 1, . . . ,j, and have 

If (yO, Pi, . . . , Pj) is in Aj or in the complement of the closure of A^, we conclude from Prop. 12.11 that 
^"■^^"^^ Tin ^°SlPSiv ((/'i.Aiv' • • • > PjAn) G 0) > -Ip,pjipi, ■■■, Pj)- 
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Thus, taking on the right-hand side the supremum over all such (pi, . . . , pj), we obtain 

The last equality uses Lemma 12.121 for the case O n Cj{p) 7^ 0, and (j2.23p is proved in this case. If 
O and Cj{p) are disjoint, then info Ii3,p,j = 00, and (|2.23p is trivially true. This ends the proof of 
Proposition [ 



2.5. The finish — proof of the LDP for {pAj^)NeN- The proof of Theorem 11.11 follows essentially 
from Proposition 12.21 and the Dawson-Gartner theorem, the LDP for projective limits, see JDZ981 
Theorem 4.6.1]. More precisely, let 

Il3,p{{pk)keN} = /3(^f{f3,p,{pk)keN) - fi/3,p) 

with 

/(/3,/5, (pfc)fegN) := sup/j(/3,p,pi, . . .,pj). 

Consider first //j^p as a function from [0, 00) to MU{oo} and endow [0, 00) with the product topology. 
By the Dawson-Gartner theorem, /^ p is a good rate function and (p^^)jvgN satisfies a large deviations 
principle with rate function /^ p. 

Now for all A^, PI^^ (PAjv ^ ^p+e) = 1- Moreover, Mpj^;, is closed as a subset of ([0, 00)^ in 
the product topology. Thus by |DZ981 Lemma 4.1.5] we conclude that {pApf)Nm satisfies a large 
deviations principle also as an Mp+e-valued random variable in this topology. 

Next, one easily sees that on Mp^^ the product topology and the £^ topology coincide. It follows 
that (paat) satisfies the LDP also in this topology with the good rate function Ijj^p. 

IjS^p is convex because it is the supremum of a family of convex functions. 

Finally, if /^,p((/Ofc)fceN) is finite, then, for all j G N, we have fj{f3,p,pi, . . . , pj) < 00 and hence by 
Proposition 12. H '^'j.^i kpk < p- Letting j — )• 00 we obtain X^^Li kpk < p. This proves that {//3,p < 00} 
is contained in Mp. 

3. Approximation with an ideal mixture of clusters 

In this section, we compare the rate function /(/3, p, •) defined in (jl.9p with an ideal rate function. 
This rate function describes a uniform mixture of clusters that do not interact with each other. This 
function has a particularly simple shape, since the combinatorial complexity does not take care of the 
excluded- volume effect, i.e., different clusters do not repel each other. 

One of the crucial points is a lower estimate for the combinatorial complexity of putting a given 
number of clusters into a large box in a well separated way. For this, we need to control the free energy 
of clusters that fit into some box of a certain volume. It is relatively easy to achieve this if the radius 
of that box is of order of the cardinality of the cluster, i.e., under the sole condition Assumption 
(V). This will turn out in Section 15.11 to be sufficient for the regime in (jl.3p . i.e., for the proof of 
Theorem II. 2 [ However, in order to handle also the much more fiexible bounds in Theorem 11.81 we will 
have to use boxes with volume of order of the cluster cardinality and to make use of Assumption 11.71 

We consider the cluster partition function, which is defined, for /3 > and /c G N, by 
Z^\l3) = ^ f e-^^^(°'^2,...,xfe)j^r|o^^2,...,Xfc} connected) dxa •• -dx^. 
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Recall the cluster partition function Z^'°(/3) with restriction to [0,a] and additional factor a~ in- 
troduced in (I2.30J) above. The reader easily checks that 

lim zf'^iP) = Zt{fi), /c G N,/3 G (0,oo). 

a— >-oo 

We also define associated cluster free energies per particle: 

ftiP) ■■= -^logZl\p), ft^\P) := -J^logZf '^(/3). (3.44) 

Let 

/^(/3) := hminf /,^l(/5) and C{P,p) := limsup/f ^^(/3), (3.45) 

fc— >oo fc->oo 

where L^ is such that the volume of [0, L^] is equal to k/p. We will see in Section [U see Lemma 14.31 
and (|4.54p . that these quantities are finite. One can actually show that they exist as limits, but we 
will not need that. 

Now we can state our bounds. The first one expresses the (simple) bound that comes from dropping 
the excluded- volume effect. Recall the definition (jl.lip of the ideal free energy f^'^'^^^. 

Lemma 3.1 (Lower bound). For all /3,p > and p G Mp, 

f{(3,p,p)>f'''-\p,p,p). (3.46) 

Proof. Recall the definition (j2.19p of the constrained partition functions Z\{(3, N, Ni, . . . ,Niy). We 
show first that 

ZA(/3,iV,iVi,...,iV^)<nM^(fll^, (3.47) 

k=i ^' 

for ah N,Ni,...,Nn G Nq with Ylk=i^^k = N. Fix such a vector {N,Ni, . . . ,Nn). Let x = 
{xi, . . . ,xn) G A^ with A^i clusters of size 1, A'^2 clusters of size 2, etc. Consider the graph with 
vertices {1, . . . , N} and edges those {i,j},i ^ j, where |xj — Xj\ < R. The graph splits into connected 
components; this induces a partition I{x) of the index set {1, . . . , A'^}. The set partition has A'^i sets of 
size 1, A''2 sets of size 2, etc. Let J = J[{Nk)k) be the collection of such set partitions of {1, ... , N}. 
Note that the integral of e~^^^ over {x : T{x) = X} does not depend on X ^ J . The cardinality of J 
is 

Therefore, for any X*"' G J ^ we may write 



ZA(/3,iV,iVi,...,Ar;v) = ^E / e-^^'-(^)l{X(rr)=X}da. 



1 



■ / e-^^^('^)]l{X(cc)=X('nd£C. 



The indicator function in the last integral can be upper bounded by dropping the requirement that 
clusters have mutual distance > R. This leads to a product of indicator functions, one for each cluster, 
encoding that the cluster is connected and stays inside A. Noting that 

J_ f e-'^^(^i'-'^'=)]l{{xi,...,Xfc} connected) dxi-.-dsfc < lAlZ^V/?) 
k\ Jan 

(integrate first over X2, . . . ,Xk at fixed xi, and then over xi), we deduce Eq. (j3.47p . 
Next, we note that n! > (n/e)" for all n G N. Therefore, ()3.47p gives that 

ZA(/3,iV,iVi,...,iV;v)<exp(-/3|A|fd-i(/3,^,(j^)^^j^)), (3.48) 
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where we have set A''^ = for A; > A^ + 1, and f^'^'^^^ is defined in (jl.lip . 

Now we turn to a lower bound for the rate function. Let O C Mp be an open set. For N G N, let 
p'^^ be a cluster size distribution in Mp of the form p^^' = A^fc/| A^vl with integer N^, and minimising 
P'^'^^^{/3,N/\A\, p) among distributions of this type. Summing Eq. ()3.48p over partitions related to O, 
we obtain 

-inf/;3,p < liminf J-logPJ^i^(pA^ G O) < -/31iminf f '^-'(/S, ^, p(-)) + /3/(/3, p). 

We have used that the number of integer partitions of A^, by the Hardy- Ramanuj an formula, is of order 
exp(0(viV)) and therefore does not contribute at the exponential scale considered here. Since Mp is 
compact, we may assume, up to choosing subsequences, that p^'^^ — ;• p^ for all k, i.e., p'^' converges 
to some p G Mp. Since the functional {p, p) i— ^ f^'^^'^^{f3, p, p) is lower semi-continuous, it follows that, 
along the chosen subsequence, 

r^-'(/3, p, p) = hminf /''i-i(/3, r^, p'"0 > ^H^ /''"'(/?, />, •)• 

We deduce 

inf/(/3,p,.)>inffd^'^i(/3,p,.), 
O o 

for every open set O C Mp. To conclude, for p G Mp, noting that Mp is metrizable, we can choose open 

environments O \ {p} and complete the proof by exploiting the lower semi-continuity of / (/3, p, •). 

D 

Our second bound controls the error when dropping the excluded- volume effect. This was much 
easier in [CKMSl^ and was hidden in the proof of Proposition 2.2 there. 

Proposition 3.2 (Upper bound). For each A; G N, let ak > be such that (a^ + R) < k/p. Then, 
for any p = {pk)k€N, 

f{l3,p,p)<^kpkfl'''''{f^)+[p-Y,kpk)r^{l3,p) + ^^Pklogp 

^ km 

Proof. We first remark that it is enough to show (j3.49p for p replaced by f e^*^ for any /c G N (where 
eC-') = [Skj)j<=N) and for p replaced by 0, the sequence consisting of zeros. Indeed, recall from 
Theorem ll.ll that /(/3, p, •) is convex, note that an arbitrary p can be written as the convex combination 



teo^-E^f-'-' + Ci-E^)". 



p V ^-^ p J 

ken ^ fceN ^ 

and note that the right-hand side of (|3.49p is afhne in p. Hence, we only have to show that 

/(/3,P,fe(^)) <p/f'^'^(/3) + i^logp+i^(-log(l-f(afc+i?)'^)+log((l+f )'^)), fc G N, (3.50) 

and that 

f{f3,p,0)<pf'^{f3,p). (3.51) 

We now prove ()3.5ip . Let O C Mp be an open set containing 0, and O its closure. By the LDP, 

hmsup -^ logPjf|^(pA^ G O) < - inf /;3,p. 

For A^ G N, consider the cluster size distribution obtained by putting all particles into one large cluster: 
pf^ = ■■■ = p'-^l^ = 0, /3^> = 1/1 A^ I . Note that p'^^ = {p''^^)km lies in Mp for any iV G N. We have 
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p(N) _s. Q g^g ^ _^ QQ g^j^(^ ^}j^jg p(N) ^ Q (2 O for sufficiently large N. As a consequence, we can lower 
bound 

^l3,AN^f^^N ^ ^) - ^i3,An^P^n - P >- Za (B N) ■ 
Recalling that \Ai\f\ = N/p, it follows that 

1 I A I 7^^^^N { fD\ 

-mU^^p > hmsup-^log ' f f l'^> > -/3prJ^{(3, p) + Pfi^, p). 

O N^oo \An\ ^Aiv(P,jVj 

Since Ip,p{-) = (3f{l3,p, •) — (3f{fi,p), this implies inf0-/(/3, p, •) < pf^{P,p). This holds for all open 
sets O containing 0. Letting O \ {0} and using the lower semi-continuity of f{(3,p,-), we deduce 

Now let us turn to (|3.5U|) . We proceed in a way analogous to Lemma 12.61 Fix fc G N. Let N be 
a multiple of k. Consider the cluster size distribution obtained by putting all particles into clusters 
of size k, i.e., put N^^^ = {N/k)5j^k for j G N. We divide the box A^r into (.^ boxes of side length 
Cfc with mutual distance at least R. Hence, i^ ~ '7;{(^k + R)~'^- The assumption (a^ + R)'^ < k/p 
guarantees that l^- > N/k for sufficiently large N. Therefore, we can lower bound 



Za,(/3, AT, Ar(-', ...,<-))> (^J^) (afZ^^"^^/?) 
Therefore, using that lA^rl = N/p and Stirling's formula. 



N/k 



liminf-^ log Za^(/3, N, N[''\ . . . , A^^') 

>^log^. (/3)-^fogP+^fog((^^^7jp-— j- 



(3.52) 



Multiplying the right-hand side with —f3~^, the right-hand side of (j3.50p arises. In the same way as 
in the proof of (j3.5ip . one derives, with the help of Lemma 12.81 that f{f3,p, f e*''') is not larger than 
— /3~^ times the left-hand side of (j3.52p . This ends the proof of (|3.50p . D 



4. Bounds for the cluster free energy 

In this section we give some more bounds that will later be used in the proofs of Theorems 11.21 and 
II. 8[ We further estimate some entropy terms, and we give bounds that control the replacement of 
temperature-depending terms by the corresponding ground-state terms. Throughout this section we 
assume that the pair potential v satisfies Assumption (V). 

We will later replace the term ^i,Pk{}ogPk - 1) in f'^''''\f3,p, {pk)k) by YjkPk^^&Pk- To this aim 
the following will be useful. 

Lemma 4.1 (Entropy bound). For any probability distribution {pk)k£N on N, 

< - ^ Pfc logPfc < 1 + log ^ kpk- 

fcGN ken 

Proof. We may assume that the expectation XlfceN ^Pk is finite. It is elementary to see that the 
maximizer of the entropy among the set of probability distributions with a given finite expectation 
is a geometric distribution. For pf^ = {I — u)u^~^ , /c e N, for some u G (0,1), the expectation is 
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X^fceN ^Pk ~ ^/(^ ~ ^) ^^^ ^^^ entropy is 

-^Pk^ogpk = -log(l -u) - {I -u)^u''~^{k- l)log?x 
feeN fceN 

U log U JT-^ u log u 



log(l --u) 



log ^ fcpfc + 

fcGN 



D 



^ Pfc log — > -2p. 



1 — u ^-^ u — 1 

ken 

We conclude by observing that xlogx > x — 1 for all x > and recalling that u < 1. 

Lemma 4.2. For any p G (0, oo) and any p = {pk)keN £ Mp, 

Pk_ 

fceN ^ 

Proof. Put m := YlkmPk ^^^^ pfc := /Ofc/"i- Then 

>^/5fclog— = y^mpk^og = mlog V m y pkiogpk 

feGN ^ fceN ' ^ keN 

m \--^ m 

> m log m — m log > A;pfc > 2m log m, 

where we applied Lemma WA\ and that X^^gN kpk < p/m. Now use the inequality xlogx > x — 1 and 
drop the term m. D 

In our bounds in Lemma [3. II and Proposition 13.21 we will later replace the cluster free energies with 
ground state energies; in this section we give bounds that will allow us to control the replacement 
error. We also prove that /^(/3) and /,^(/3, p) are finite. 

Lemma 4.3 (Lower bound for f^{P) and /^(/3)). There is a constant C > such that for all 
/3g(0,oo), 



In 



C 



ken,(3 e (0,oo). 



particular, f^{l3) > Coo — % for any (3 G (0,oo). 



Proof. We follow [ClKMSlOl Sec. 2.4]. First, note that 



1 



^,^k/3)<e-'^^'=- {(x2,...,x,)E 



i,d\fc-l 



A;! 



: {0, X2, ...,Xfc} i?-connected} 



with I • I the Lebesgue volume. Now, with each x' = (x2, . . . , Xk) such that x := (0, x') is -R-connected, 
we can associate a tree T{x') with vertex set {1, . . . , fc} and edge set E{T{x')) C {{i, j} : i 7^ j}, and 
such that 

{i,j} G ^(r(a:')) ^ \xi -Xj\< R. 

Note that for a given a;', there are in general several trees satisfying this condition; we pick arbitrarily 
one of them and call it T(x'). Now we have 

{x' G (M'^)'=-^ I {0,x') /^-connected} 
= X] k^' ^ (M'^)^"^ I {0,x') i?-connected, T{x') = T} 

T tree 

- Z] |{^' ^ {M.'^f-^ I (0,cc') i?-connected,{i,i} G E{T) 



Xi\ < i?} 



r tree 
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For each given tree T, the Lebesgue volume of the set in the last line above is upper bounded by 
\B{0,R)\^~^. By Cayley's theorem, see |AZ981 pp. 141-146], the number of labeled trees with k 
vertices is A;'^"^. Thus 

and the proof is easily concluded. D 

Now we show that the volume constraint in the cluster partition function is immaterial for large /? 
if the radius of the confining box is of order of the particle number with a sufficiently large prefactor. 

Lemma 4.4 (Low-temperature behavior of f^'°'{P)). For any fc G N and any choice of ak{f3) in 
[kR, oo), 



lim/f"'=('^)(/3) = :|. 

Proof. The lower bound, '>', is trivial since Z^'"(/3) < Z^{P) for any o. For afc(/3) > kR, the box 
[0,afc(/3)]'^ is certainly large enough to contain a minimiser of £C i— >• Uk{x). Therefore, lower bounding 
the integral by an integral in a neighborhood of the minimiser, we find 

hmmf-logZ, ' ''^'^' > — — , 
/3^oo /3 " k 

which is the upper bound '<'. D 

Under additional assumptions, most importantly Assumption II. 7^ it will be enough to pick a^ of 
order k^'"^ instead of k, with some error of order 4 log/3: 

Lemma 4.5 (Uniform low-temperature bounds for f^'^{f3)). Suppose that the pair potential also 
satisfies As sumptions \1.6\ and \1.7\ There is an a > and a (3 > such that for all j3 £ [13, oo), and 
every sequence of au's satisfying a^ > ak^''^, 

/f"'^-(/3)<^ + ^log/5, ken. (4.53) 

In particular, for any p G (0, 1/a'') and (3 G [/3,oo), 

ft{P,p)<eoo + ^\ogf3. (4.54) 

Proof. The strategy of the proof is as follows. According to Assumption 11.71 'we may pick a minimiser 
for Uk that fits into some ball whose volume is of order of the particle number. Then we restrict the 
integral in the definition of the cluster partition function to some neighbourhood of this minimiser 
and control the error with the help of the Holder continuity from Assumption 11.61 Let us turn to the 
details. 

Let c > be as in Assumption 11.71 (^ > as in Lemma 12.51 Then a := 2(c + 5) satisfies ak^''^ > 
6 + ck^''^ for all /c G N. Fix t G {l,R/b). Let nmax G N be the maximal number of particles that can 
be placed in B{0,R), keeping mutual distance > r^am, with rmm as in Assumption 11.61 

For A; G N, let a^ > ak^''^ and let x^°^ = {xi , . . . ,x^^) be a minimiser of the energy C/^ that fits 
into the cube with side length a^ — 5. Thus x^"^ is 6-connected, and \xi — Xj\ > r^i^ for every i / j. 
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The scaled state tx'-"^ is t6-connected and has minimum interparticle distance > trmin- By the Holder 
continuity of the potential v, 

k 

\U{tx^^^) - U{x^"^)\ <\Y^ Y\v{t\xf' - xf\) - v{\xf' - xf 

i=l j^i 

< A;nmaxsup{|v(r') - v{r)\ : r > r^am, r' > rmm, \r - r'\ < {t - 1)6} 

with C and s such that \v{r') — v{r)\ < C\r' — r\'^ for any r,r' > rmin- Let e G (0, 1) such that 

e < (5/2, Tmin < trmin — 2e, and tb + 2e < R. 

We will obtain a lower bound for Z"^ '""^ (/3) by considering configurations (xi, . . . , x^) with exactly one 
particle per e-ball around tx^"^ for j = 2, . . . ,k. To this end, put 

M':= U {B{tx^:i2ys)x...xB{tx^:i,^,e)] 

where &'i^_i denotes the set of permutations of 2, . . . ,k, and let Ai be the set of configurations in 
the cube of side length ak — 6 obtained by rigid shifts from configurations in {a;^° } x A4' . For small 
enough e, the balls B{tx'-^Ls , e) , . . . , B{tx'^}^. , e) do not overlap, and M.' has therefore Lebesgue volume 

(A:- l)!|S(0,e)|^-^ Moreover, 

\M\ > \M'\ {at- 5- ck^l'^f > ^\M'\. 
Now a; G A^ is i?-connected and has minimum interparticle distance > rmin- Thus 

\U{x)-U{tx^°^)\<Ckn^^^e\ x e M. 
Restricting the integral in the definition (j2.30p of Z^ '"'°(/3) to A^, we obtain 

4zf''''={f3) > |||S(0,e)|'=-iexp(-/3(£;fc + CA:n^axe^)). 



This implies, for |S(0,e)| < 1 

ft'^HP) < ^ + ^'"T^ - ^log \B{0,e)\ + 



ci,a,(^^ < _^ + J- - -log \B{0,e)\ + -^. 



Now we pick e = 1//3 for definiteness and obtain that (I4.53P is satisfied for sufficiently large /3. D 

5. Proof of F-convergence and uniform bounds 

In this section, we prove Theorems 11.21 and 11.81 Recall that Theorem 11.21 is proved under the sole 
Assumption (V) and that we additionally suppose that Assumptions 1 1 . 6l and 11.7] hold for Theorem 1 1.8i 

5.1. Proof of Theorem 11.21 Fix v G (0, oo) and let (0, oo) 9 s i— )■ (/3(s), p{s)) be a curve in (0, oo)^ 
such that, as s — )• oo, 

/3(s)^oo, p{s)^0, --— -log p{s) ^u. 

We need to show that, for any q = (gfc)fcGN G Q, 
Lower bound: For all curves q'°^ — t- q, 

liminf ^/(/3(s),p(s),p(^)) > g,{q). (5.55) 
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Upper bound / recovery sequence: there is a curve q'"' — )• q such that 

limsup-^/(/3(s),p(s),p'=0 < g,{q). (5.56) 

Ptoo/ 0/ i/ie /ower bound. We write q''' = {ql^^)k S Q- Define p^"' = {p^j^^)km by g^"^ = kp'j^^/p. Let 
C > and^> such that kff{l3) > Ek-Ck(3~^ for any k G NU{oo} and /3 G [^,00), see Lemma IJTSl 
Then Lemma |3. II gives 



^^ le 






((-..^ Z^ nC*^ V ^ 



/3(s)f^P(^)^ P(^) ^ /3(s 



The term in the second line converges to gvisi) because of the continuity of the map q 1— )■ X^fcgpj QkiEk — 
v)/k + (1 — 'Ylk&'i1k)^oo] here enters the property E^/k — >■ 600- The terms in the last line are, by 
Lemma l4.2t of order l//3(s) and therefore converge to 0. D 

Proof of upper bound / existence of a recovery sequence. We choose p-dependent box sizes ak{p) such 
that {ak{p) + RY < k/{2p), a^ > R, and a^ > 6 + k^''^{r]^c + ^), with 6 as in Lemma [231 Such a choice 
is possible for small enough p, and compatible with the additional requirement that ak{p) — )■ 00 as 
p — )• 0, for every /c G N. Lemma 12.51 tells us that 

which can be upper bounded by some constant C, uniformly in A; G N and sufficiently large s. 
Now we apply Prop. [3?2l This gives, for sufficiently large s and any sequence p = {pk)k, 

^fWisiPishP) < E k^fi'-^^^^^^^ Wis)) + (1 - E ^^) r^i^^pi^)) 

Consider first the case ^1^=1 Ik ~ ^- ^^^ 1^"^ '■~ 9- ^^ have, for any K £ N, 

'^^ ' k=l '^ ' k=K+l ^^ ' 

Since ak{p{s)) — )• 00 as s — )• 00 for any fc G {1, . . . , K}, using Lemma [^^ we get 

limsup^^/(/3(5),/,(s),p'=))< Vgfc^4-^ + C V Qk- 

'^^ P^'^' fc=l '' k=K+l 

Letting K ^ 00 we find that limsnpg^^^ p{s)^^ f{l3{s), p{s), p'^''^) < g,y{q). 

Next, consider the case qk = for all A; G N. For n G N, let s^ > large enough so that for s > Sn, 
|/n'"" — En/n\ < 1/n. The sequence {sn)nm can be chosen increasing and diverging. We set 

k{s) := n for s G [s„, Sn+i) and n G N. It follows that k{s) — )• 00 as s — )■ 00, and 



30 SABINE JANSEN, WOLFGANG KONIG, BERND METZGER 

from which we deduce /^A-)*" (/3(s)) — )• eoo as s — )• oo. Set qk{s) := 5^ ^(s). Then we find 

hmsup— -/(/3(s),p(s),p<"') < Coo = gu{q)- 

s— 7>oo P{S) 

To conclude, we observe that every q G Q can be written as a convex combination of a vector q' 
with X^fcgp^ 9fc = 1 ^'^d the zero vector, and a recovery sequence is constructed by taking the convex 
combination of q' and the recovery sequence for the zero vector. D 

5.2. Proof of Theorem 11.81 Proof of (1): We prove (jl.lSp in terms of p^'s instead of qts. Then it 
reads 

f{l^,pAPk)k&i) -Y^Pk\Ek + —^j + [p-^kpkje^ <— plog/3, (pfc)fceN G Mp. (5.58) 

Lemmas 13. H 14.21 and 14.31 vield that there is C € (0, oo) such that, for all (3,p £ (0,oo) and p = 

f{(3,p,p)>f'''''\P,p,p) 

fcGN ^ fcGN '^ ^ fceN '^ ^ fcGN 

fceN ^ fcGN " 

(5.59) 
This is '>' in (j5.58p . For proving '<', we pick, for p G (0,cx3) and /c G N, box diameters ak{p) such 
that ak{p) > ak^' , with a as in Lemma 14.51 ^ind {ak{p) + i?) < k/2p, for all /c G N. This is possible 
provided p < k/2{ak^' + R) for any fe G N, and this is, by monotonicity in k, guaranteed for p <p, 
where we put p = 2(0+ bW • We may also assume, without loss of generality, that a > R, which implies 

that afc(p) > R for all /c G N. We obtain, for {/3,p) G [/3,oo) x (0,p), and C > as in Lemma [i3| for 
any p G Mp, with the help of Proposition 13.21 

f{P,P,p)<Y.kpJ^-^logf^) + (p-Y,kp,)(e^-^logp) + ^PY.P^ 



fcgN fc6N ^ ^ fceN 



+ ^2^'°'^(-^°s(i-^) + i°^ 



< Y^Pk^Ek + i^) + (p - 5Z ^/Ofc)eoo + ^log/3 + (t^+ 1)^ log 2, 

which is the corresponding upper bound in (15.580 . 

Proof of (2): Let p = {pk)k be a minimiser of f{(3,p, •) and q := {kpk/ p)k£n- Write u = — /3~^ log p. 
Then 

-/(/3, P) = -/(/?, /5, P) > 9u{q) - § log /? > p(i^) - § log /?. 
Similarly, let g be a minimiser of gu{') and p := {pqk/k)keN- Then 

Pii^) = 9u{q) > -fW, P^P)-J log l3>-f{f3,p)-j log /?. 

Proof of (3): Let p = (pfc)fc be a minimiser of f{/3,p, •) and g := {kpk/p)keN- Write 1^ = —/3~^ log p. 
Then (1) and (2) yield 

gAq) - ^H < -fW,P,p) + ^iog/3 - (-/(/?,p) - ^iog/5) < 2^iog/5- 
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Hence, 

2^ log /? > gM) - ^i{u) = Y, {^^ - K^))Qk + (eoo - K^)) (l - J] Qk) • (5.61) 

For 1/ < I/*, we use that /i(i^) = eoo and estimate 

Ek- V _ r ^ _ Ek - fcCoo - V ^ V* - L' 

rC fv fv 

Substituting this in (I5.6ip . this yields the first claim, (ll.lTp . 

For 1/ > I/* , we restrict the first sum on the right of (j5.6ip to k £ N\M(i/), where we lower estimate 
the brackets against A(z^), and we estimate eoo — /^*(^) ^ ^(^)- This gives 

2|log/3> Yl A(i/)g, + AH(l- J]gfc)=A(z.) Y Ik- 

^ k<m\M{u) ken keM(u) 

This yields the second claim, (jl.lSp . 

6. Appendix: Proof of Lemma 11.31 

Here we prove Lemma 11.31 With the exception of the positivity of ly*, this has been proved in 
[CKMSl^ Theorem 1.5]; that proof works under the slightly different assumption on v that we have 
here. To obtain the positivity of u*, this proof needs a slight modification, which we briefly indicate 
now. Fix M,N £N. Let tc'^' = (xi, . . . , xn) £ (M'^)^ be a minimiser of Un and y'-"'^ = (yi, . . . , vm) 
a minimiser of Um- Recall that b is the potential range and let 5 > be such that v < on (b — 6,b). 
Let e e (0, (5/2). Let a G M'^ be such that the shift y^"^ := (yi, . . . , y^/) '■= ivi + O', ■ ■ ■ i Vm + cl) satisfies 

• all points from y'*^' and ic'^' have distance \xi — yj\ >b — 5 + e (and hence v{\xi — yj\) < 0), 

• there is at least one pair of particles {xi,yj) with distance \xi — yj\ <b — e. 

Let ;e(^+^^' := (a;(~),y<''>) G {R'^)^+^^ . Let c := - sup^.^[i^__s+s,b-e] ^(^) > 0- Then we have 

E^+M < Uix^'^+^y) < C/(a;(^)) + U{y^''') -c = En + Em - c. 
In particular, the sequences {Ej\f)j\f^fq and {Ej\f — c)]\f^^ are subadditive, whence 

T En .. En - c En — c 
eoo = hm — — = lim — — = mt — — . 

Af^oo N N^oc N NeN N 

Because of the stability of the pair potential, we have Coo > — oo. The inequality eoo < {En — c)/N 
for any A^ leads to En — Ne^o ^ c for any A, and this is the positivity of v* . 

Acknowledgement: We gratefully acknowledge financial support by the DFG-Forschergruppe 718 
'Analysis and Stochastics in Complex Physical Systems.' 



References 

[AZ98] M. AlGNER and G.M. Ziegler, Proofs from THE BOOK. Springer, Berlin (1998). 

[CKMSIO] A. COLLEVECCHIO, W. KONIG, p. MORTERS and N. SiDOROVA, Phase transitions for dilute particle systems 
with Lennard- Jones potential, C'ommun. Math. Phys. 299:3, 603-630 (2010). 

[CLY89] J.G. CONLON, E.H. Lieb and H.-T. Yau, The Coulomb gas at low temperature and low density, Commun. 
Math. Phys. 125, 153-180 (1989). 

[DZ98] A. Dembo and O. Zeitouni, Large Deviations Teehniques and Applications. 2nd edition, Springer, Berlin 

(1998). 



32 



SABINE JANSEN, WOLFGANG KONIG, BERND METZGER 



[DS84] R. DiCKMAN and W.C. Schieve, Proof of the existence of the cluster free energy, J. Stat. Phys. 36, 435-446 

(1984). 

[F85] C.L. Fefferman, The atomic and molecular nature of matter. Rev. Mat. Iberoamencana 1, 1-44 (1985). 

[GHMOl] H.-O. Georgii, O. Haggstrom and C. Maes, The random geometry of equilibrium phases, Phase transi- 
tions and critical phenomena. Vol. 18, Academic Press, San Diego, CA, 2001, pp. 1-142. 

[H56] T.L. Hill, Statistical Mechanics: Principles and Selected Applications, McGraw-Hill Book Co., Inc., New 

York (1956). 

[HLOl] J.-B. Hiriart-Urruty and C. Lemarechal, Fundamentals of Convex Analysis, Springer, Berlin (2001). 

[LP77] J.L. Lebowitz and O. Penrose, Cluster and percolation inequalities for lattice systems with interactions, 

J. Stat. Phys. 16:4, 321-337 (1977). 

[dM93] G. Dal Maso, An Introduction to T -convergence, Birkhauser, Boston (1993). 

[M75] M.G. MiJRMANN, Equilibrium distributions of physical clusters, Commun. Math. Phys. 45:3, 233-246 (1975). 

[PY09] E. Pechersky and A. Yambartsev, Percolation properties of the non-ideal gas, J. Stat. Phys. 137:3, 

501-520 (2009). 

[R81] C. Radin, The ground state for soft disks, J. Stat. Phys. 26 (1981), 365-373. 

[R99] D. Ruelle, Statistical Mechanics: Rigorous Results. World Scientific, Singapore (1999). 

[S03] N. Sator, Clusters in simple fluids, Phys. Rep. 376, 1-39 (2003). 

[T06] F. Theil, a proof of crystallization in two dimensions, Comm. Math. Phys. 262, 209-236 (2006). 

[YFS09] Y.A. Yeung, G. Friesecke and B. Schmidt, Minimizing atomic configurations of short range pair poten- 
tials in two dimensions: crystallization in the WulfF shape, preprint, arxiv:0909.0927vl (2010). 

[Z08] H. Zessin, a theorem of Michael Miirmann revisited, Izv. Nats. Akad. Nauk Armenii Mat. 43:1, 69-80 

(2008), MR 2465000 (2010c:60149). 



